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ABSTRACT It is shown that Schroedinger equation is not consistent with information theory. From the 
modified form of information which ensures that the most probable density function it yields tallies with a 
general form of continuous Riemann integrable density function that has real or imaginary zeros or 
singularities at end points of [a,b]eR, a new variational formulation for quantum mechanics is proposed 
that yields a system of Euler-Lagrange equations that are non-linear. It is proved that the solutions of this 
system are unique, orthonormal and complete. One dimensional harmonic oscillator has been solved. 

I INTRODUCTION 

Quantum mechanics makes probabilistic predictions and information theory (Shannon et 
al, 1969 and Jones, 1979) yields the most probable density function subject to a given set 
of constraints (Venkataraman, 1999). If the probabilistic predictions of quantum 
mechanics (Messiah, 1968 and Schiff, 1968) are also the most probable ones, then 
Schroedinger equation will be consistent with information theory. An affirmative answer 
to the title question would then render a comparative study of an ensemble of classical 
systems, generated, for instance, by replacing the initial conditions of classical equations 
with random variables, with quantum system meaningful. In section II the question is 
answered in the negative and a new variational formulation of quantum mechanics is 
proposed in section III. The latter involves modified form of information 
(Venkataraman, 1999) that ensures the probability density function of a real variable 
obtained from it tallies with a general form of continuous Riemann integrable density 
function in [a,b] that has real or imaginary zeros or singularities at end points. In 
appendix 1 multivariate probability densities and, convergence of power series in several 
variables using spherical polar coordinates are discussed. In appendix 2 it is proved that 
for the proposed variational formulation the Lagrange multipliers corresponding to 
orthogonality are zero and that the orthonormal set of solutions is unique and complete in 
the corresponding Hilbert space. In section IV the system of non-linear Euler-Lagrange 
equations is solved for one dimensional harmonic oscillator. Results are summarised in 
section V. 

II. Information Theory, Probability and Schroedinger Equation 

Information theory predicts the most probable density function subject to a given set of 
constraints. If p(x)dx is the probability density in [a,b]czR, the information / associated 
with it is given 



where the brackets denote the averaging done with respect to p(x)dx. If the constraints 
are 
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where c;'s are constants, then p(x)dx is obtained by minimizing / subject to the given 
linearly independent set of constraints. The variational method yields 



p(x)dx = cxp(-Lf^ ax l )dx 

where di's are Lagrange multipliers to be determined from 



(3) 



= Ci 



(4) 



Schroedinger equation can also be derived from the variational principle of minimizing 
(H) where H is the Hamiltonian subject to normalisation of eigenvectors. (lnp(x)) does 
not occur in the variational formulation and the inference thereof is that Schroedinger 
equation is not consistent with information theory or that the predictions of quantum 
mechanics are not the most probable ones. 

Before providing an alternative variational formulation of quantum 
mechanics, a natural question to ask is "Is the probability density function specified in 
equation (3) by fixing all moments the most general one in [a,b] ? " If p(x)dx is 
continuous and has no real or imaginary zeros or singularities in [a,b] Taylor expansion 
of Inp (x) about xo^[a,b] yields 

In p(x) = In p(xo) + £°° (x-Xof 1 / i ![^ln p(x)] (x . n) (5) 

dx l 

From the above equation, equation (3) follows, as its R.H.S. is convergent. Hence 
equation (3) represents a continuous probability density function that has no zeros or 
singularities in [a,b]. If p(x)dx has real or imaginary zeros or singularities at end points 
of [a,b] specified by Z p and S p , R.H.S. of equation (5) will not be convergent but 
Taylor expansion of ln{p/(Z p S p )} will converge yielding 

p(x)dx = ZpS P Qxp(-i:jL ax l )dx (6) 

If p(x)dx has finite discontinuities or real or imaginary or zeros or singularities in [a,b] at 
a finite number of interior points then [a,b] has to be subdivided into a finite number of 
intervals in each of which an expression of form (6) will hold. This form of density is 
non-negative in the entire domain unlike the Gram-Charlier series that seeks an 
improvement on the Gaussian. It also explains why central limit theorem yields a 
Gaussian. Equation (6) will ensure faster convergence of approximate solutions than a 
power series when differential equations are solved (Venkataraman, 1999). In the case of 
more than one variable there will be more constraint equations specified for example in 
two variable case by 

(x l y J ^ = Cij (7) 

which yields 

p(x, y)dxdy=exp[-Z~j )=(0 ,o) a^xy 3 ]dxdy (8) 

where {x,y)sDczR and p does not vanish anywhere in D. Coefficients a,y are derived 
from 



d 

{InpO)} ) = dj (9) 



Vide the discussion in appendix 1. Probability densities of the form (6) for one variable 
could be derived from information theory if the latter is redefined as 

I = (\n{p/(Z P S P )} ) (io) 



III . Quantum Mechanical Equation Consistent with Information Theory : 

Since information given by equation (10) yields a probability density that allows zeros 
and singularities at end points of [a,b] and that p(x)dx being given by \yA 2 dx in quantum 
mechanics where symbol \ff is the state vector, zeros of f can also be in the interior of 
the domain, the following variational formulation is proposed: 
Minimize / given by (10) subject to (H) being constant or vice versa which yields 

- 1 / 2d 2 / dx 2 + V (x) = l[l + In { W *¥„ I (ZTZ) } ]¥„ 

(ii) 

in the place of Schroedinger equation, where Z„ represents the zeros of \j/„ in D, the 
domain of integration and \f/„eL 2 (D) and A« is the Lagrange multiplier. These equations 
are nonlinear. The logarithmic term, depending upon the sign of the Lagrange 
multiplier, will increase/decrease the energy values. If normalisation is also introduced 
in equation (11) Schroedinger equation can be retrieved in the limit of the Lagrange 
multipltier corresponding to the logarithmic term tending to zero. The solutions of (11) 
are the normalised solutions. Since the above problem is the same as minimising energy 
subject to fixing the information, energy values will also be ordered as in the standard 
linear case. 

In appendix 2 it is shown that the multipliers corresponding to orthogonality are zero and 
that equation (11) yields the unique, complete set of solutions in the appropriate Hilbert 
space. 

The reader who recognizes the importance of using all moments as constraints 
in deriving equation (3) using (2) and (1) might ask why (ft) for i>l have not been 
concluded in the variational formulation . At present the linear equation 

[H+pH*T¥ = M (12) 

arising from specifying (H 2 ) is being studied. Inclusion of the latter will yield 
correction terms. There are papers in literature where the logarithmic no n- linear terms 
are introduced and the equations studied with different aims (Doebner et al, 1999, 
Weinberg, 1989). 

IV One Dimensional Harmonic Oscillator: 

Equations (11) reduce to 

- 1 / 2 d 2X ¥ n I dx 1 + x 2 1 2^ n = L[l + In / (Z*Z) } ] v P n 

(13) 

where the zeros of y/„ are represented by a polynomial P„ of degree n. Starting from 



x i f n (x) = x k P n (x)exp[-Y,2 a2jc ] (14) 

and taking recourse to the differential equation it can be shown that 

= K(Jljix) exp [-a, - fix 2 ] (i5a) 



k=0orl (15b) 

A. = (4jSL 2 -l)/4j& (15c) 

/3L 2 = [2a - 1] / [8(n + + OL)] (i5d) 

exp(2o,) = VTl2 M ft!/ (16) 



where //„'s are Hermite polynomials . The solutions to (16) being unique y/n's given by 
(15) are unique. One can alternatively prove that the solutions are complete using the 

2 2 

standard basis in L ® : If y/eL ® then it can be written as a linear combination of the 
vectors of the standard basis, the set of Hermite functions It is easy to show then the 
orthogonality of this with those given in (15a) implies it is a null vector. The energy 
values in the various states y/„ are given by 

K = Ajl- 2a, - (2n + 1) / 2] (i?) 

The numerical values of various quantities are given below for two cases: 



N 






K 


E„ 





0.561903 


0.165957 


-1.34046 


0.836186 


1 


0.8846183 


0.182575 


-1.18673 


2.69296 


2 


1.483947 


0.265717 


-0.675132 


3.01642 


3 


2.374767 


0.271151 


-0.650844 


4.71831 


4 


3.3791495 


0.312319 


-0.488143 


5.00752 


5 


4.5328009 


0.309387 


-0.498664 


6.76468 


6 


5.7558755 


0.334322 


-0.413460 


7.03368 


7 


7.07846158 


0.330258 


-0.426725 


8.81483 



Since A/'s are negative the energy values are higher than in the linear case though the 
spacing between energies is not constant.. 

V Conclusion : A modified form of information is proposed in equation (10) so that the 
probability density obtained from it tallies with a continuous Riemann integrable one 
that could have real or imaginary zeros or singularities at end points of [a,b]. For the 
case of many variables the method of information theory can be used to derive joint 
probability density as in equations (7) and (8). Power series expansion, equivalent to 
Taylor series, for two variables is given in appendix for standard functions. To ensure 
consistency of quantum mechanics with most probable predictions, the new variational 
formulation (35), that yields a unique, complete orthonormal set of solutions is suggested 
. As the proof in appendix 2 shows to obtain the orthonormal set of solutions in the 
Hilbert space to nonlinear Schroedinger equation one has to take recourse to the 
variational method to obtain the other conditional minima.. The non- linear Euler- 
Lagrange equations have been solved analytically for one dimensional harmonic 
oscillator. The results could be compared with experimental results only when the 



problem is solved in L (-1,1) as in reality the systems are always finite (Venkataraman, 

1999) preferably after incorporation of (H 1 ) for at least i=2 in (1 1). 
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Appendix 1 

A probability density of the form given in equation (8) can be used when solutions to 
differential equations are sought for. Before deriving it in yet another way power series 
in several real variables is considered below. 

It is well known (Courant et al 1989, Vol.1) that starting from 

Pn(x) = Y,l aJc l (18) 

where 

a = 1/ illd 1 / dx' P(x)](* = o) (19) 

Taylor and his pupils defined power series starting from 

/ M = Lio * H\d l I dx l f{x)\, = 0) + Rn (20) 

In the above equation it has to be shown that the remainder Rn tends to zero to justify 
taking the upper limit to infinity. 

Analogously if one begins from 

P^y) = HH M) a iJ x i y j (2i) 

and defines 

an - -J— [<9< i+ ''> / dx'dy^l, = <>,, = <» (22) 

one could show, since 

(1+ axf = 1+ k(ax) 1 1!+ k(k - l)(ax) 2 1 2!+ (23) 

is convergent for lax I < 1 for a and k real 



(1+ X y) k = 1+ k(xy) 1 1!+ k(k - l)(xy) 2 1 2!+ (24) 

is convergent for \xy\<l or equivalently r<7 in polar coordinates, (x,y)eR 2 . Similarly 
since the power series for exp(a x) is convergent for all xeR. where a is real 

exp(xy) = l+(xy)/l!+(xy) 2 /2!+ +(xy) r I >!+.... (25) 

and the above series is convergent for all {x,y)eR . Similarly the power series 
expression for logarithm could also be derived to arrive at equation (8). Taylor 
expansion for two variables is defined by 

The coefficients in (24) and (25) and (26) are precisely those derived from an equation 
analogous to (22). Hence power series in two variables can be defined through 

where 

I 

an = [d {i+j) / dx l dy j f] ( x = <>,, = <» (28) 

It is easier to consider convergence questions for functions in equations (8), (24), (25) 
and (28) of several variables using polar coordinates. Equation (27) is the limiting form 
of (26) and gives a more natural expression for ay 's. Maxima and minima could also be 
found easily solving just one equation in the radial coordinate and checking the signature 
of the second derivative in that coordinate (R.P.Venkataraman 1999). 
APPENDIX 2 
The equation 

m=Wo (29) 
corresponds to finding the lowest eigenvalue of H subject to the normalisation 

V¥o=l (30) 



\dP¥ 



in which D is the domain of integration and dl is the infinitesimal volume 
element. The equation for y/i arising out of the additional constraint of 
orthogonalisation with respect to y/o is 

tf¥i = Aii¥i + AJo¥o (31) 
where jlo is the Lagrange multiplier corresponding to the condition of 
orthogonality. It can be shown that /n is zero using self adjointness of H and 
the orthogonality condition. Thus the eigenvalue Xi is the next conditional 
minimum. For any n it can be shown that 

= Ln^n (32) 

because the Lagrange multipliers corresponding to the orthogonality with ty/yt 
for all k<n are zero and that the eigenvalues are ordered. The purpose of the 
above is to show that the solution of variational problems are unique and that 
equation (32) is a special case in which the system of equations gets reduced 
to one. 



Starting with a variational formulation consistent with information theory 
the system of nonlinear Euler-Lagrange equations 

- 1 / 2 d 2X ¥ J dx 2 + V(xf¥. = l[l + In {M, / (Z*Z) } pP, 

(33) 

where A„ is the Lagrange multiplier and Z„ represents the zeros of was derived 
by minimising 

/ = ( ln{¥,*¥, /(Z*Z)} ) (34) 

subject to being constant or vice versa. For any variation 8\[/„* in \ffn 

jdT(W + #F.*)[-ln{ + &¥n*y¥» I (Z*Z) } + (l/L)Hf¥ n 

D 

- jdr^:[-\n{^*^n/(Z*Z n )} + (l/l,)Hfi' n = (35) 

D 

which yields equation (33) for y/„.. It shall now be proved that the Lagrange 
multipliers corresponding to orthogonality are zero. Choosing n=l one obtains 
introducing orthogonality with \f/o 

- 1 / 2 d 2X ¥ J dx 2 + V(xy¥n = Ml + lnjW^ ,/ (Z*Z) } pP, + 

(36) 

Multiplying the above equation by y/o and integrating , one obtains for the multiplier 

\iX \ ar\dvm^:^« = j jt[ w ]Aoin{ (^0*^0) / (z *z ) 



- JrfnWi^Po]Ailn{(^i^i)/(Z*Z)}^iW (37) 
D 

where all the terms have been multiplied by the factor in square brackets to facilitate 
integration by parts. Using \[/o8\[/ =S( \f/ 2 /2) \[/i8\f/i =5 ( \f/i 2 /2) one obtains on using the 
finiteness of the logarithmic term and the boundary conditions of the solutions the value 
of 

Ho to be zero. To prove the multipliers corresponding to orthogonality for all \jf„ are 
zero one assumes the result true for k=0 to (n-2) and proves as above that the multiplier 
corresponding to k=(n-l) is zero. If the nonlinearity is of the form ly/2f then it can be 
shown easily the multipliers corresponding to orthogonality are nonzero. 

To prove that the solutions of (35) are unique, let the contrary be assumed. Let 
y/i and be two normalised solutions corresponding to the given value of A. Then the 
energy values of these two vectors could be found. It is easy to show as above that 

Er(X)-E 2 (X) = (38) 

which is a contradiction. 



As mentioned in the text, this formulation gives the solutions corresponding to the 
complete set of conditional minima of the energy functional and hence 

EoK Eu..K Ei.. To prove that this orthonormal set is complete let the contrary be 
assumed. Let y/eL (D) be a nonzero vector and let {(ft} a basis in that space such that 
{(ffVjeCfiD) . Then yf can be written as a linear combination of the latter and the 
expectation value of the Hamiltonian operator in this state can be found. Then this value 
E must lie between two values of the energies, say Ei and E(i+i). This is a contradiction 
as E(i+i) is the next conditional minimum after Ei. Hence the solutions corresponding to 
the complete set of minima of the energy functional must be complete in the 
corresponding Hilbert space. 



